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Abstract 

We propose a method for the treatment of two-point boundary 
value problems given by nonlinear ordinary differential equations. The 
approach leads to sequences of roots of Hankel determinants that con- 
verge rapidly towards the unknown parameter of the problem. We 
treat several problems of physical interest: the field equation deter- 
mining the vortex profile in a Ginzburg-Landau effective theory, the 
fixed-point equation for Wilson's exact renormalization group, a suit- 
ably modified Wegner-Houghton's fixed point equation in the local 
potential approximation, a Riccati equation, and the Thomas-Fermi 
equation. We consider two models where the approach does not apply 
in order to show the limitations of our Pade-Hankel approach. 
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1 Introduction 



Some time ago Fernandez et al [1-9] developed a method for the accurate 
calculation of eigenfunctions and eigenvalues for bound states and resonances 
of the Schrodinger equation. This approach is based on the Taylor expansion 
of a regularized logarithmic derivative of the eigenfunction. The physical 
eigenvalue is given by a sequence of roots of Hankel determinants constructed 
from the coefficients of that series. One merit of this approach, called Riccati- 
Pade, is the great convergence rate in most cases and that the same equation 
applies to bound states and resonances. Besides, in some cases it yields upper 
and lower bounds to the eigenvalues [1]. 

The logarithmic derivative satisfies a Riccati equation, and one may won- 
der if the method applies to other nonlinear ordinary differential equations. 
The purpose of this paper is to investigate whether a kind of Pade-Hankel 
method may be useful for two-point boundary value problems given by non- 
linear ordinary differential equations. 

In Section [2] we outline the method, in Section [3] we apply it to several 
problems of physical interest, and in Section H] we discuss the relative merits 
of the approach. 



It is our purpose to propose a method for the treatment of two-point bound- 
ary value problems. We suppose that the solution f{x) of a nonlinear ordi- 
nary differential equation can be expanded as 



about X = 0, where a and (3 are real numbers, and (3 > 0. We also assume 
that we can calculate sufficient coefficients fj in terms of one of them that 
should be determined by the boundary condition at the other point; for 
example, at infinity. We show several illustrative examples in the following 
section. 

We try a rational approximation to x^"f{x) of the form 
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where z = x" . The Taylor expansion of the usual Pade approximant yields 
M + N + 1 coefficients of the series ([T]) [11]; but in the present case we require 
that the rational approximation ([2]) gives us one more coefficient, that is to 
say, M + + 2. If M = + ci, = 1, 2, . . ., d = 0, 1, . . ., this requirement 
leads to the equation [1-9] 

Hi, = \ fi+j+d+i\ij^Q^i^_^ = 0, (3) 

where D = + 1 = 2, 3, . . . is the dimension of the Hankel determinant Hf). 

In general, equation ([3]) exhibits many roots and one expects to find a 
sequence, for D = 2, 3, . . . and fixed d, that converges towards the required 
value of the unknown coefficient. From now on we call it Hankel sequence for 
short. If such convergent sequence is monotonously increasing or decreasing 
we assume that it yields a lower or upper bound, respectively. Such bounds 
where proved rigorously for some eigenvalue problems [1]. 



3 Examples 

In order to test the performance of the Pade-Hankel method, in this section 
we consider the examples treated recently by Boisseau et al by means of a 
most interesting algebraic approach [10]. We first consider the field equation 
determining the vortex profile in a Ginzburg-Landau effective theory [10] 
(and references therein) 

fir) + i/'(r) + (l - ^) fir) - f{rf = 0, r > 0. (4) 

The solution /(r) satisfies the expansion ([1]) with x = r, a = n = 1,2, . . ., 
and (3 = 2. If we substitute this series into the differential equation and solve 
for the coefficients fj, we obtain them in terms of the only unknown /o that 
is determined by the boundary condition at infinity: f{r ^ oo) = 1 [10] (and 
references therein). The coefficients fj, and therefore the Hankel determinant 
Hi, are polynomial functions of /q. For example, for n = 1 we have 

f = -i^ f = Jr— f- = ■^Q _ ^li. (^) 
■'^ 8 ' 192 24' ■''^ 9216 576''" ^ ^ 

Tables [T] and [2] show two Hankel sequences with d = and d = 1 that 
converge rapidly towards the result of the accurate shooting method [10] for 
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n = 1 and n = 2, respectively. We appreciate that in the case n = 1 the 
sequences with d = and d = 1 give upper and lower bounds, respectively, 
that tightly bracket the exact value of the unknown parameter of the theory: 
0.58318949586060 < /o < 0.58318949586061. 

On the other hand, the appropriate Hankel sequences are oscillatory when 
n > 2 and their rate of convergence decreases with n. Table [3] shows the best 
estimates of /o for n = 2,3, 4. 

Our second example is the fixed-point equation for Wilson's exact renor- 
malization group [10] (and references therein) 

2/" (x) - 4/(x)/'(x) - 5x/'(x) + fix) = 0, X > 0. (6) 

The solution to this equation can be expanded as in equation ([1]) with a = 1 
and (3 = 2. The first coefficients are 

For large values of x the physical solution should behave as f{x) = ax^^^ + 
a^/(5x'^/^) + . . .. The Hankel sequences with d = and d = 1 converge 
towards the numerical result [10] (and references therein) from above and 
below, respectively. Fig. [1] displays the great rate of convergence of these 
sequences as A = \fo{D, d = 0) — fo{D, d = D = 2,3, . . ., from which we 
obtain the accurate bounds -1.22859820243702192438 < fo < -1.22859820243702192437 

The third example comes from a suitably modified Wegner-Houghton's 
fixed point equation in the local potential approximation [10] (and references 
therein) 

2r(x) + [1 + f (x)][5/(x) - xf'ix)] =0, x>0. (8) 

The solution satisfies the series ([T]) with a = 1 and (3 = 2, and the first 
coefficients are 

3 3 ' ~ 60 15 60 ^ ' 

On the other hand, the acceptable solution should behave as f{x) = ax^ — 
4/(3x) + . . . when x ^ 1. 

Table H] shows Hankel sequences with d = and d = 1 that clearly 
converge towards the numerical value of /o [10] (and references therein). 

We have also applied our approach to the ordinary differential equation 
for the spherically symmetric skyrmion field [10] (and references therein) but 
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we could not obtain convergent Hankel sequences. We do not know yet the 
reason for the failure of the method in this case. 

Present approach has earlier proved suitable for the treatment of the 
Riccati equation derived from the Schrodinger equation [1-9]. Consider, for 
example, the following Riccati equation 

f'{x)- f{xf + x^ = 0, x>0. (10) 

The solution can be expanded as in equation ([1]) with a = P = 1; the first 
coefficients are ^ 

/l = /o ) /2 = /o5 /s = ~2 + /o ) • ■ ■ 

There is a critical value /oc of /(O) = /o such that f{x) ~ —x at large x 
if /(O) < /oc, f{x) develops a singular point if /(O) > /oc, and /(x) ~ x 
at large x if /(O) = /oc- Present Pade-Hankel method yields the value of 
/oc with remarkable accuracy as shown in Table O The rate of convergence 
of the Hankel sequence for this problem is considerably greater than for the 
preceding ones. 

If we substitute f{x) = —y\x)/y{x) into equation ffTOl) . then the function 
y{x) satisfies the Schrodinger equation for a harmonic oscillator with zero 
energy on the half line: y"{x) — x'^y{x) = 0, and the problem solved above is 
equivalent to finding the logarithmic derivative at origin y'{0)/y{0) so that 
y{x) behaves as exp(— at infinity. Obviously, any approach for linear 
differential equations is suitable for this problem. 

Finally, we consider the Thomas-Fermi equation [11-13] (and references 
therein) 

$"(a;) = x~'^^^{xf/\ x>0 (11) 

that provides a semiclassical description of the charge density in atoms of 
high atomic number. It poses the problem of finding the slope at origin 
$'(0) so that $(0) = 1 and $(x oo) = 0. The change of variables t = x^^'^ 
leads to a more tractable equation 

te'{t)-Cit)-4t'm'^',t>0 (12) 

where ^(t) = ^(t^). We can expand the solution to this equation in a Taylor 
series: ^(t) = 1 + a2t^ + 4t^/3 + 2a2t^/5 + . . . and the unknown slope is given 
by the unknown coefficient: $'(0) = a2- One easily derives a recurrence 
relation for the coefficients of this expansion [13]. Although the treatment 
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of equation fll2p is straightforward, we find it more convenient to define the 
function f{t) = ^(t)^/^ that is a solution to 

fit)nt) + f{tf] - fit)f\t) - 2t'fitf = 0. (13) 

The function f(t) satisfies a series like ([1]) with a = and (3 = 1. The first 
coefficients are 

/l=0, /3 = ^, /4 = -Y, /5 = -^,... (14) 

We expect that a Hankel sequence will converge towards /2 = /"(0)/2 = 
^'{0)/2. Since /4 is the first nonzero coefficient that depends on /2 we choose 
d > 4 in the Hankel determinant. Table [6] shows that the Hankel sequence 
with d = 4 converges rapidly giving a most accurate value of the slope at 
origin $'(0). 

Finally, we consider two examples discussed by Bender et al [14]; the first 
of them is the instanton equation 

f"{x)+f{x)-f{xf = (15) 

with the boundary conditions /(O) = 0, /(oo) = 1. The solution to this 
equation is f{x) = tanh (^x/y/2j. The expansion of f{x) is a particular case 
of equation ([1]) with a = 1 and /3 = 2; its first coefficients being 

. ^_fo fo + 1) . /o (66/o^ + 1) 

Ji g , 120 ' 5040 ' ■ ■ ■ ' ^ ' 

where fo = /'(O) is the unknown. The Hankel series with d = and d = 1 
converge rapidly giving upper and lower bounds, respectively, to the exact 
result /o = l/\/2. 

The second example is the well known Blasius equation [14] 

2y"'{x)+y{x)y"{x) = (17) 

with the boundary conditions y{0) = y'{0) = 0, y'(oo) = 1. The expansion 
of the solution in a Taylor series about x = is a particular case of equation 
([1]) with a = 2 and /3 = 3; its first coefficients are 

f = -— f = ^^-^Q fl8l 
•'^ 60' 20160'' ^ ^ 

Since, in general, fj oc /q^^, then the only root of the Hankel determinants 
is /o = that leads to the trivial solution y{x) = 0. We thus see another 
case where the Pade-Hankel method does not apply. 
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4 Conclusions 



We have presented a simple method for the treatment of two-point boundary 
value problems. If there is a suitable series for the solution about one point, 
we construct a Hankel matrix with the expansion coefficients and obtain the 
physical value of the undetermined coefficient from the roots of a sequence 
of determinants. The value of this coefficient given by a convergent Hankel 
sequence is exactly the one that produces the correct asymptotic behaviour 
at the other point. We cannot prove this assumption rigorously, but it seems 
that if there is a convergent sequence, it yields the correct answer. More- 
over, in some cases the Hankel sequences produce upper and lower bounds 
bracketing the exact result tightly. 

Present Pade-Hankel approach is not as general as the one proposed 
by Boisseau et al [10], as we have already seen that the former does not 
apparently apply to the skyrmion problem or to the Blasius equation [14]. 
However, our procedure is much simpler and more straightforward and may 
be a suitable alternative for the treatment of this kind of problems. Besides, 
if our approach converges, it yields remarkably accurate results as shown in 
the examples above. For example, it gives us the slope of the electrostatic 
potential of the Thomas-Fermi theory with unprecedented accuracy. 



References 






[1] 


Fernandez 


F 


M, Ma Q, and Tipping R H 1989 Phys. Rev. A 39 


[2] 


Fernandez 


F 


M 


1992 Phys. Lett. A 166 173. 


[3] 


Fernandez 


F 


M 


and Guardiola R 1993 J. Phys. A 26 7169. 


[4] 


Fernandez 


F 


M 


1995 J. Phys. A 28 4043. 


[5] 


Fernandez 


F 


M 


1995 J. Chem. Phys. 103 6581. 


[6] 


Fernandez 


F 


M 


1995 Phys. Lett. A 20 3 275. 


[7] 


Fernandez 


F 


M 


1996 J. Phys. A 29 3167. 


[8] 


Fernandez 


F 


M 


1996 Phys. Rev. A 54 1206. 


[9] 


Fernandez 


F 


M 


1997 Chem. Phys. Lett 281 337. 



7 



[10] Boisseau B, Forgacs P, and Giacomini H 2007 J. Phys. A 40 F215. 

[11] Bender C M and Orszag S A 1978 Advanced mathematical methods for 
scientists and engineers (McGraw-Hill, New York). 

[12] Plindov G I and Pogrebnya S K 1987 J. Phys. B 20 L547. 

[13] Fernandez F M and Ogilvie J F 1990 Phys. Rev. A 42 149. 

[14] Bender C M, Pelster A, and Weissbach F 2002 J. Math. Phys. 43 4202. 



8 



Table 1: Convergence of the Hankel series for the connection parameters of 
the global vortex for n = 1 



D 


d 


= 


d 


= 1 


2 


0, 


.595 


0. 


,578 


3 


0, 


.584 


0. 


,5829 


4 


0, 


.58324 


0, 


,58315 


5 


0, 


.58320 


0. 


,583183 


6 


0, 


.583192 


0. 


,583187 


7 


0, 


.583190 


0. 


,5831890 


8 


0, 


.5831897 


0. 


,5831893 


9 


0, 


.58318954 


0. 


,58318946 


10 


0, 


.58318952 


0. 


,58318948 


11 


0, 


.58318951 


0. 


,583189491 


12 


0, 


.583189498 


0. 


,583189494 


13 


0, 


.5831894964 


0. 


,5831894953 


14 


0, 


.5831894961 


0. 


,5831894956 


15 


0, 


.5831894960 


0. 


,5831894957 


16 


0, 


.58318949590 


0. 


,58318949583 


17 


0, 


.58318949588 


0. 


,58318949584 


18 


0, 


.583189495867 


0. 


,583189495854 


19 


0, 


.583189495864 


0. 


,583189495857 


20 


0, 


.583189495862 


0. 


,5831894958591 


21 


0, 


.5831894958609 


0. 


,5831894958598 


22 


0, 


.5831894958607 


0. 


,5831894958601 
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Table 2: Convergence of the Hankel series for the connection parameters of 
the global vortex for n = 2 



U 


a 


n 

= U 


a 


1 

= i 


o 
o 


0. 


,iOO 


0. 


,ioi 


4 


0. 




0. 


,io4 


r 




0. 


iOcSiO 


0. 


,iooO 


D 


U. 


iooUy 


A 
U. 


,iooii 


/ 


U. 




0. 


1 c:o 1 n 
loolO 


Q 
O 


U. 


,it)oUyy i 


U. 


,it)oiU 


n 

y 


U. 


.ioouyyi 


U. 


,ibouyy 


10 


0. 


,15309914 


0. 


,1530989 


11 


0. 


,15309912 


0. 


,153099095 


12 


0. 


,15309917 


0. 


,153099091 


13 


0, 


,153099105 


0, 


,153099097 


14 


0. 


,1530991021 


0. 


,15309911 


15 


0. 


,15309910272 


0. 


,153099102 


16 


0. 


,153099102697 


0. 


,153099103 


17 


0. 


,153099102782 


0. 


,15309910292 


18 


0. 


,153099103124 


0. 


,15309910293 


19 


0. 


,153099102857 


0. 


,15309910289 


20 


0. 


,153099102864 


0. 


,15309910278 


21 


0. 


,15309910286136 


0. 


,153099102860 


22 


0. 


,15309910286142 


0. 


,153099102858 



Table 3: Best estimates of the connection parameters of the global vortex for 
n = 2, 3, 4 by means of Hankel sequences with D < D^^x 



n- Dmax fo 

2 21 0.15309910286 

3 21 0.0261834207 

4 26 0.0033271734 
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Table 4: Convergence of the Hankel sequences for the Wegner-Houghton 
connection parameter 



D 




(i = 




d= 1 


3 


-0, 


.3013652092 


-0, 


.4190129312 


4 


-0, 


.5405112824 


-0, 


.4696457170 


5 


-0, 


.4552012493 


-0, 


.4604796926 


6 


-0, 


.4624525979 


-0, 


.4616935821 


7 


-0, 


.4613759926 


-0, 


.4615091717 


8 


-0, 


.4615571129 


-0, 


.4615373393 


9 


-0, 


.4615303767 


-0, 


.4615331535 


10 


-0, 


.4615342975 


-0, 


.4615338165 


11 


-0, 


.4615336147 


-0, 


.4615337043 


12 


-0, 


.4615337357 


-0, 


.4615337227 


13 


-0, 


.4615337173 


-0, 


.4615337196 


14 


-0, 


.4615337207 


-0, 


.4615337202 


15 


-0, 


.4615337200 


-0, 


.4615337201 


16 


-0, 


.46153372013 


-0, 


.461533720119 


17 


-0, 


.461533720113 


-0, 


.4615337201157 


18 


-0, 


.4615337201168 


-0, 


.4615337201163 


19 


-0, 


.4615337201161 


-0, 


.4615337201162 


20 


-0, 


.4615337201162 
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Table 5: Convergence of the Hankel sequences with d — ior the Riccati 
equation. 



D fo 

4 0.6762 

5 0.675970 

6 0.6759785 

7 0.67597823 

8 0.6759782403 

9 0.675978240059 

10 0.6759782400675 

11 0.675978240067277 

12 0.6759782400672850 

13 0.675978240067284722 

14 0.675978240067284729 

15 0.67597824006728472899 

16 0.67597824006728472900 

17 0.67597824006728472900 




Figure 1: A = \fo{D, d — 0) — fo{D, d—l)\ for Wilson's renormalization 
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Tabic 6: Convergence of the Hankel sequences with d = 4 for the Thomas- 
Fermi equation 



D 




2/2 


10 


-1 


.5880709 


11 


-1 


.5880706 


12 


-1 


.58807103 


13 


-1 


.588071024 


14 


-1 


.5880710227 


15 


-1 


.58807102264 


16 


-1 


.588071022609 


17 


-1 


.588071022609 


18 


-1 


.5880710226116 


19 


-1 


.5880710226115 


20 


-1 


.58807102261139 


21 


-1 


.58807102261138 


22 


-1 


.58807102261137 


23 


-1 


.58807102261137 


24 


-1 


.5880710226113756 


25 


-1 


.58807102261137537 


26 


-1 


.58807102261137532 


27 


-1 


.5880710226113753154 


28 


-1, 


.5880710226113753152 


29 


-1 


.5880710226113753154 


30 


-1 


.5880710226113753137 
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